Evidence for nodal quasiparticles in electron-doped cuprates from penetration depth 

measurements 
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The in-plane magnetic penetration depth, A(T), was measured down to 0.4 K in single crystals of 
electron-doped superconductors, Pri.ssCeo.isCuCU-i (PCCO) and Ndi.ssCeo.isCuO^a (NCCO). 
In PCCO, the superfluid density varies as T 2 from 0.025 up to roughly 0.3 T/T c suggestive of a 
d-wave state with impurities. In NCCO, A(T) shows a pronounced upturn for T < 4 K due to the 
paramagnetic contribution of Nd 3+ ions. Fits to an s-wave order parameter over the standard BCS 
range (T/T c = 0.32) limit any gap to less than A mi „(0)/T c = 0.57 in NCCO. For PCCO, the absence 
of paramagnetism permits a lower temperature fit and yields an upper limit of A m in{0)/T c — 0.2. 



There is by now a consensus that the hole-doped high- 
T c cuprates exhibit wave pairing symmetry [fl| [l0| . 
For electron-doped cuprates [0 the issue remains un- 
resolved. While most theories for the mechanism of high 
temperature superconductivity are insensitive to the sign 
of the carriers, some predict that n and p type materials 
will have different pairing symmetry, making its determi- 
nation an important challenge |l^,|l3|. Early microwave 
measurements of the penetration depth in NCCO were 
interpreted within an s— wave model p4-16|. However, 
Cooper pointed out that the power law dependence for 
A(T) indicative of a nodal order parameter could be 
masked by a large paramagnetic contribution from Nd +3 
ions Jirj l . Newer microwave measurements by Kokales 
et ai, performed on the same sample used in this pa- 
per, have revealed an upturn and power-law tempera- 
ture dependence and are consistent with our data Jl8| . 
Measurements of A(T) using single grain boundary junc- 
tions [jl9| have favored a gapped state. Some tunneling 
measurements favor an s— wave order parameter, albeit 
with significant departures from an isotropic weak cou- 
pling BCS picture |^o| , while others report a zero-bias 
conductance peak ]21-23|. Half-integral flux indicative 
of d-wave pairing pi was recently reported in tricrystal 
experiments with both NCCO and PCCO films @. 

In this letter we report measurements of A(T) down to 
0.4 K in single crystals of both NCCO and PCCO. Lower 
temperatures and higher resolution combine to permit a 
more precise determination of the temperature depen- 
dence of A(T) than any previously reported. In NCCO, 
a large paramagnetic contribution is observed below 4 K. 
In non-magnetic PCCO, we find an overall T 2 variation 
of the superfluid density up to T/T c w 0.3, suggesting the 
presence of nodal quasiparticles in the presence of strong 
impurity scattering. 

Single crystals of Ri.85Ceo.isCu04-5 (R=Nd or Pr) 



were grown using directional solidification technique and 
annealed in argon to achieve optimal superconducting 
properties [^5| . Penetration depth was measured using 
an 11 MHz tunnel-diode driven LC resonator [^6] psf . 
Samples were mounted on a movable sapphire stage with 
temperature controllable from 0.4 K to 100 K. The low 
noise level, Af min /fo «5x 10~ 10 , results in a sensitivity 
of AA < 0.5 A for our samples [0.5 x 0.5 x 0.02 mm]. 
The large anisotropy of these materials (X c /X a b ~ 30-80 
]l5] | ) forces one to apply the rf field perpendicular to the 
conducting planes. Otherwise, the frequency shift will 
be dominated by changes of the interplane penetration 
depth, for which there exists no straightforward connec- 
tion to the pairing symmetry. A semi-analytical solution 
for the rf susceptibility of a platelet sample of square base 
2w x 2w and thickness 2d in this orientation was analyzed 
in detail in Ref. p8[ . At low temperatures the frequency 
shift, A/ = f{T) — /(0), is related to the change in pen- 
etration depth, AA = A(T) - A(0), via A/ = — GAA, 
with the calibration constant G = V s fo/ [2Vb(l — N)R], 
where N is the effective demagnetization factor, V s is 
the sample volume, Vq is the effective coil volume, /o 
is the resonance frequency in the absence of a sample 
and R w 0.2w is an effective dimension |28 . Although 



this result is similar to the known solution for an infinite 
slab in parallel field || , the effective dimension R differs 
significantly from R = w/2 obtained for an infinite bar. 
This difference is due to penetration of the magnetic field 
from the top and bottom surfaces. The sample and ap- 
paratus dependent constant A/ = V s fo/ [2Vq(1 — N)] 
is measured by moving the sample out of the coil in situ. 
The overall calibration was tested with samples of Nb, 
YBCO and BSCCO and gave dX/dT within 10 % of 
reported values |sj . In order to determine the normal- 
ized superfluid density, p s = [A(0)/A(T)] 2 , it is neces- 
sary to know the absolute magnitude of the penetration 
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depth, A(0). Measured values of A(0) in NCCO vary from 
1000 - 2600 A 0l||2|g. In PCCO the only reported 
value, A(0) w lOOOA, was estimated from the measure- 
ments of the lower critical field H c \ and is less reliable 
due to demagnetization and possible vortex effects fl3l| . 
We recently developed a new technique for determining 
A(0) from the frequency shift obtained by warming a 
sample coated with a thin layer of Al above T c {Al). A 
detailed description of this procedure will be published 
elsewhere |32 . This technique applied to PCCO gave 
A(0) = 2500 ± lOOl which will be used as an upper 
estimate in this paper. 
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FIG. 1. Low temperature variation of the penetration 
depth AA(T) in PCCO single crystal. Lines show fits to dif- 
ferent models described in the text. 

In Fig.|l| we plot the penetration depth for PCCO 
sample 1 along with several fits. The fitting range 
was 5 K to assure validity of the low temperature 
BCS expansion for an isotropic s-wave state, AA = 
A(0) t/ttA(0) /2T exp (- A(0) /T) Q. In each case, the 
small negative offset A — A(0) — A(0.4 K) was deter- 
mined as a fit parameter. The solid line is a fit with 
188 data points to a power law, AA = A + BT 2 with 
B = 3.70 ± 0.01 A/K 2 and X 2 ~ 8.2. The short dot- 
ted line shows the best fit to the BCS s-wave expres- 
sion. With both A(0) and A(0) as free parameters we 
obtained A(0)/T c = 0.538 ± 0.002 and A(0) = 426 ± 3 
A. The s-wave fit is somewhat worse than the power law 
(x 2 ~ 17.6) and gives an unrealistic value for A(0). The 
dash-dotted line shows the s-wave fit where A(0)/T c was 
fixed at the weak-coupling BCS value (1.76). In this case 
an unrealistically large A(0) = 13570±50 A was obtained. 
For comparison, the dotted line is a plot of the BCS ex- 
pression with A(0)/T c = 1.76 and set to a more realistic 
value of A(0) = 2500 A. 

If the order parameter is an anisotropic s-wave, then 
the minimum gap value determines the low temperature 
asymptotic behavior. The BCS functional form for AA 
still holds, but with A m i„ replacing the isotropic gap. 
The temperature range over which this asymptotic form 



is valid is now reduced accordingly. For an isotropic 
gap, A(0)/T c = 1.76, the range of validity in reduced 
temperature is (T/T c ) max = t max w 0.32 0. For an 
anisotropic gap, simple rescaling forces the range of va- 
lidity down to t max ~ 0. 18 A m i„(0)/T c . Without a pri- 
ori knowledge of A m i„(0)/T c , we do not know t max and 
so it is necessary to successively reduce the range until 
the gap value obtained from the fit becomes independent 
of the range. Following this procedure for PCCO, with 
A(0) = 2500 A fixed, we find that A(min) /T c extrapo- 
lates to 0.20 ± 0.05 as t max — > 0. The same procedure 
for Nb yields A(0)/T c = 1.74 ± 0.02, as expected. This 
stricter criterion would imply that any residual gap is less 
than 11% of the isotropic BCS value. 

The overall best fit for sample 1 was achieved for a 
y2.25±o.oi p 0wer i aw (^2 ^ 2 Although this may ap- 
pear unphysical, it is important to recall that the integer 
power laws expected for a nodal order parameter strictly 
apply to the normalized superfluid density, p s , and not 
the measured quantity, AA. If p s = 1 — c n T n then AA 
has corrections of order T 2n and a fit to AA can result in 
an artificial intermediate power (e.g. 2.25). This distinc- 
tion is clearly evident in high quality, untwinncd YBCO 
above 10 K |27|l. 
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FIG. 2. Superfluid density p a vs. (T/T c ) 2 , calculated as- 
suming A(0) = 1000, 1550 and 2000 A. Dotted lines indicate 
pure quadratic power law. Dash-dotted line shows isotropic 
s— wave behavior. 

In Fig. H we plot, for sample 1, normalized superfluid 
density, p s w [1 + AA(T)/A(0)p 2 , vs. (T/T c ) 2 for three 
choices of A(0) spanning the range of reported values. 
A(0) =1200 A yields 1 -p s ~ (T/T c ) 2 up to 8.4 K, while 
larger choices for A(0) reduce the range of pure quadratic 
behavior. For comparison we show data taken for Nb in 
the same apparatus. Up to to T/T c — 0.5, the Nb data 
fits the low temperature BCS expansion perfectly with 
A(0)/Tc = 1.74 ± 0.02, giving us confidence that the 
measurement technique is sound. 

With an exponent of n = 2.25, sample 1 is our weak- 
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est candidate for a nodal order parameter. Fig. |3| shows 
data for samples 1, 2 and 3 (offset for clarity) and the 
same Nb data with a fit to the BCS form. Samples 2 
and 3 have power laws much closer to n = 2. In this plot 
we have chosen the largest value of A(0) in order to cast 
the power law model in the most unfavorable light (i.e., 
smallest range of pure quadratic behavior). We conclude 
that the superfluid density in PCCO is best described by 
a quadratic power law variation with temperature. 
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FIG. 3. p a in three PCCO single crystals assuming 
A(0) = 2500 A (data offset for clarity). Dotted lines show 
T 2 fits at low temperatures. Also shown data for Nb and 
calculated isotropic s-wave curve. 

There are currently two theories for a quadratic power 
law in d-wave superconductors. Kosztin and Leggett 
showed that the divergence of the effective coherence 
length near the nodes of a d— wave order parameter yields 
1 — p s ~ T 2 due to nonlocal electrodynamics Q . Non- 
locality is predicted to arise for the orientation used 
here (field perpendicular to conducting planes) below 
Tnoniocai ~ £(0) A(0)/A(0), where £(0) is the coherence 
length. In electron-doped cuprates T non i oca i « 0.5 K to 
2.5 K within our current knowledge of superconducting 
parameters. Since we observe a quadratic temperature 
dependence up to 8-10 K in some samples, nonlocality 
is unlikely to be the source. A stringent test for nonlo- 
cality would require a comparison between this data and 
A(T) obtained from the H\\ab plane orientation. How- 
ever, the H\\ab orientation involves the interplane pene- 
tration depth, as discussed earlier. 

Impurity scattering in the unitary limit provides a 
more plausible explanation for the quadratic dependence 
of p s (T) |||]. In the "dirty d-wave" scenario, p s (T) will 
cross over from a linear to quadratic temperature depen- 
dence below T* ~ 61n27/7r, where 7 ~ 0.63^TA(0) 
and r is a scattering rate parameter, proportional to the 
impurity concentration The slope d\/dT 2 \ T ^ Q ~ 



7rA(0)/ (67A(0)). Casting this result in dimensionless 
form gives: T/T c ~ 0.28/ [~dp s /d{T/T c ) 2 ] . In Fig. | 
we plot T c versus T/T c for all five samples studied. Sam- 
ples 1,2 and 3 are marked. The transitions in some sam- 
ples were broad and two different criteria were used to 
estimate T c - onset of the diamagnetic signal and the in- 
flection point of the AA(T) curve. In general, the trend 
shows T c suppression with increased scattering rate as 
expected for a d-wave state with impurity scattering [|| . 
r/T c is at least 10 times larger than the scattering rate 
observed in clean YBCO §J7|. 



26 




18 1 ' 1 ' 1 ' 1 

0.1 0.2 0.3 0.4 

T/T 

c 

FIG. 4. T c versus scattering parameter F/T c determined 
using 5 PCCO crystals. Filled symbols show T° nset denned 
as the onset of diamagnetism. Open symbols show T™ de- 
fined as the inflection point on AA(T). 

We now discuss measurements of A(T) in NCCO. This 
compound has been studied much more thoroughly than 
PCCO and was cited as the first evidence for s— wave 
pairing in e-doped materials [ pTp5[ |. However, Nd 3+ ions 
introduce a large paramagnetic background and influence 
the measured penetration depth Jj"7| . With magnetic per- 
meability pi(T) = 1 + C/ (9 + T), the measured pene- 
tration depth is given by A(T) = Al(T) y/ /i(T), where 
Al is the London penetration depth, C is a Curie- Weiss 
constant, and O is the characteristic temperature for an- 
tiferromagnetic interaction JTJJ. To fit the data,we take 
= 1.2 K from neutron scattering |m| and specific heat 
p5| measurements. For C we have chosen two represen- 
tative values C = 0.3 and 0.05, calculated assuming the 
effective magnetic moment of Nd 3+ ions to be 2Aps 0] 
and \\iB-, respectively. 

Figure || shows AA(T) measured in a single crystal 
of NCCO. The inset shows the low-temperature range. 
Below T w 4 K there is a pronounced upturn, which 
we attribute to the paramagnetic contribution of Nd 3+ 
ions. The upper solid line in the inset to Fig.|| shows the 
power law fit (Al oc T n ) which yields n = 1.35 ±0.03 and 
n = 1.40 ± 0.03 for C = 0.3 and 0.05, respectively. (Fits 
for the two different values of C are indistinguishable on 
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this scale.) The lower line shows analogous fits to the low 
temperature s— wave expression from which we obtain 
A(0)/T c = 0.569 ± 0.006 and A(0)/T c = 0.573 ± 0.006. 
Fits were obtained from data up to t = 0.32. The higher 
temperature data is shown for completeness. Changing 
the value of from 1.2 to 2 also had a small effect on the 
fit parameters. The value of A(0)/T c is close to that ob- 
tained in PCCO. Again, for a strict test we should fit over 
a correspondingly reduced temperature range. However, 
the dominant paramagnetic contribution below t < 0.16 
renders this procedure meaningless. The n = 1.4 ex- 
ponent obtained from the power law fit is closer to the 
clean d-wave limit and could imply that after correction 
for paramagnetism, unitary-limit scattering in NCCO is 
smaller than in PCCO. 
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FIG. 5. AA(T) in NCCO single crystal. Inset: Low tem- 
perature range. Lines are the fits assuming a power law or an 
s— wave temperature dependence for the London penetration 
depth Xl(T) as described in the text. 

In conclusion, we have measured the penetration depth 
A(T) in electron-doped PCCO and NCCO single crys- 
tals down to 0.4 K. In non-magnetic PCCO, p s decreases 
quadratically with temperature up to t « 0.3, consistent 
with a dirty d-wave (gapless) scenario. The correlation 
between T c and the rate of change of supcrfluid density 
is also consistent with this picture. In NCCO, a large 
paramagnetic contribution to the penetration depth was 
observed. Aj,(T) was found to vary as T 1A . For both 
materials, a fit over the same temperature range to an 
s-wave model sets an upper limit of A(0)/T c = 0.57 but 
requires unrealistically small values A(0). For PCCO, the 
test can be made more stringent and we reduce the upper 
limit to A mi „(0)/T c = 0.2. 
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